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ABSTRACT 

Quadratic  forms  on  complex  random  matrices  and  their 
joint  eigenvalue  densities  are  derived  for  applications  in  in¬ 
formation  theory.  These  densities  are  represented  by  com¬ 
plex  hypergeometric  functions  of  matrix  arguments,  which 
can  be  expressed  in  terms  of  complex  zonal  polynomials. 
The  derived  densities  are  used  to  evaluate  the  most  impor¬ 
tant  information-theoretic  measures,  the  so-called  ergodic 
channel  capacity  and  capacity  versus  outage  of  multiple- 
input  multiple-output  (MIMO)  Rayleigh  distributed  wire¬ 
less  communication  channels.  Both  (spatially)  correlated 
and  uncorrelated  channels  are  considered  and  the  corre¬ 
sponding  information-theoretic  measure  formulas  are  de¬ 
rived.  It  is  shown  how  channel  correlation  degrades  the 
communication  system  capacity. 

1.  INTRODUCTION 

Let  an  n  x  m  complex  Gaussian  random  matrix  X  be  dis¬ 
tributed  as  X  ~  CA(0,Si  ®  £2)  with  mean  £{X}  =  0 
and  covariance  cov{X}  =  £1  ®  £2,  where  £1  £  <Cnxn 
and  £2  €  Cmxm  are  positive  definite  Hermitian  matrices. 
Then  the  quadratic  form  on  X  associated  with  the  positive 
definite  Hermitian  matrix  A  is  defined  by 

S  =  XH  AX. 

Here,  we  study  the  distribution  of  S,  denoted  by 
CQn,m(A,  £1,  £2),  and  its  application  to  information  the¬ 
ory.  We  also  derive  the  joint  eigenvalue  densities  of 

S,  which  are  represented  by  complex  zonal  polynomials. 
Complex  zonal  polynomials  are  symmetric  polynomials  in 
the  eigenvalues  of  a  complex  matrix,  see  [5],  and  they  en¬ 
able  us  to  represent  the  derived  densities  as  infinite  series. 
If  A  =  In,  £1  =  In  and  £2  =  £,  then  S  =  XHX 
is  said  to  have  a  complex  Wishart  distribution,  denoted  by 
CWm(n,  £),  see  [3],  [6]  and  references  therein. 


The  theory  of  quadratic  forms  on  complex  random  ma¬ 
trices  is  used  to  evaluate  the  capacity  of  multiple-input, 
multiple-output  (MIMO)  wireless  communication  systems. 
Let  us  denote  the  number  of  inputs  (or  transmitters)  and 
the  number  of  outputs  (or  receivers)  of  the  MIMO  wireless 
communication  system  by  nt  and  nr,  respectively,  and  as¬ 
sume  that  the  channel  coefficients  are  distributed  as  com¬ 
plex  Gaussian  and  correlated  at  both  the  transmitter  and  the 
receiver  ends.  Then  the  MIMO  channel  can  be  represented 
by  an  nrxnt  complex  random  matrix  H  ~  CN(0,  £r®£t), 
where  £r  and  £t  represent  the  channel  correlations  at  the 
receiver  and  transmitter  ends,  respectively.  This  means  that 
the  covariance  matrices  of  the  columns  and  rows  of  H  are 
denoted  by  £r  and  £t,  respectively.  If  £r  =  In„  (or  <r'1 2 * * S,  IHr ) 
and  £4  =  <y2Int  (or  Int)  then  the  channel  is  called  uncor¬ 
related  Rayleigh  distributed  channel.  The  information  pro¬ 
cessed  by  this  random  channel  (or  mutual  information  of 
this  random  channel)  is  a  random  quantity  which  can  be 
measured  in  two  ways,  namely,  ergodic  capacity  (or  aver¬ 
age  mutual  information)  and  capacity  versus  outage  (or  x 
percent  outage). 

Recent  studies  show  that  a  MIMO  uncorrelated 
Rayleigh  distributed  channel  achieves  ahnost  n  more  bits 
per  hertz  for  every  3dB  increase  in  signal-to-noise  ratio 
(SNR)  compared  to  a  single-input  single-output  (SISO)  sys¬ 
tem,  which  achieves  only  one  additional  bit  per  hertz  for 
every  3dB  increase  in  SNR.  where  n  =  inin{ut,  nr},  see 
[2]  and  [12],  In  this  context,  we  may  mention  a  MIMO 
strategy  that  offers  a  tremendous  potential  to  increase  the 
information  capacity  of  single  user  wireless  communica¬ 
tion  systems,  namely,  the  Bell-Labs  Layered  Space-Time 
(BLAST)  architecture,  see  [8]  and  references  therein.  But 
in  many  practical  situation,  channel  correlation  exists  due 
to  poor  scattering  conditions,  which  degrades  capacity  (see 
[1],  [9],  [10]  and  references  therein).  In  [1],  the  authors 
studied  the  correlated  channel  matrix  with  large  dimension 
(asymptotic  analysis),  which  is  only  an  approximation  to  the 


practical  correlated  channel  matrix  with  finite  dimension.  In 
[10],  the  effect  of  channel  correlations  on  MIMO  capacity  is 
quantified  by  employing  an  abstract  scattering  model.  More 
recently,  in  [9]  the  authors  have  derived  an  exact  ergodic 
capacity  expression  for  an  uncorrelated  Rayleigh  MIMO 
channels,  which  is  different  from  the  work  of  Telatar  [12], 
and  have  also  derived  an  upper  bound  to  ergodic  capacity 
for  correlated  Rayleigh  MIMO  channels.  In  this  paper,  we 
first  derive  the  densities  of  quadratic  forms  on  complex  ran¬ 
dom  matrices  and  their  joint  eigenvalue  densities.  Then, 
using  these  densities  we  evaluate  the  capacity  degradation 
for  the  correlated  channel  matrices  H  ~  CN( 0,  Sr  <g>  St) 
with  nr  >  Tit,  nr  =  nt  and  nt  >  nr.  This  will  be  done  by 
deriving  closed-form  ergodic  capacity  and  outage  capacity 
formulas  for  correlated  channels  and  then  numerical  evalu¬ 
ation.  Note  that  this  work  is  an  extension  of  our  early  work 
[6],  where  we  have  extensively  studied  the  complex  Wishart 
matrices  and  uncorrelated/correlated  (only  at  the  transmitter 
end)  MIMO  channel  capacities. 

This  paper  is  organized  as  follows.  Quadratic  forms  on 
complex  random  matrices  are  studied  in  Section  2.  The  ca¬ 
pacity  of  a  MIMO  channel  and  the  computational  methods 
are  given  in  Sections  3  and  4,  respectively. 


2.  QUADRATIC  FORMS  ON  COMPLEX  RANDOM 
MATRICES 

In  this  section,  the  densities  of  quadratic  forms  on  complex 
random  matrices  are  given  and  then  joint  eigenvalue  den¬ 
sities  are  derived.  The  probability  distributions  of  random 
matrices  are  often  derived  in  terms  of  hypergeometric  func¬ 
tions  of  matrix  arguments.  In  the  sequel,  we  need  to  use  the 
following  complex  hypergeometric  function  of  two  matrix 
arguments, 


o  F^\X,Y)  =  ^ 
k= 0  k 


Ck(X)Ck(Y) 

k\CK(In) 


(1) 


and  X[k](X)  is  the  character  of  the  representation  [k]  of  the 
linear  group  given  by 


xmPO  = 


det  [(A^+m"i)] 

[(*r')] 


as  a  symmetric  function  of  the  eigenvalues,  Ai, . . .  ,  Xm,  of 
X.  Note  that  real  and  complex  zonal  polynomials  are  partic¬ 
ular  cases  of  (general  a)  Jack  polynomials,  ( X ),  where 

a  =  1  for  complex,  and  a  =  2  for  real,  zonal  polynomials, 
respectively.  In  this  paper  we  only  consider  the  complex 
case;  therefore,  for  notational  simplicity  we  drop  the  super¬ 
script,  a,  of  Jack  polynomials,  as  was  done  in  equation  (2), 
i.e.,  we  write  CK(X)  :=  (X).  Finally,  we  have 


CK(In)  =  2 2kk!  (  i 
where 

( 

2 


!(H„ 


H.-6G 


n:<j(2ki-2kj~i+j) 


rii=i(2*j  +  r  -*)! 


,  (3) 


^(n-i  +  1)) 

/  ki 


and  the  partition  k  of  k  has  r  nonzero  parts.  Here  (a)k  = 
o(o  +  1)  •  •  •  (o  +  k  —  1)  and  (o)o  =  1. 

The  next  theorem  gives  the  density  of  quadratic  forms 
on  complex  random  matrices  S  =  XHXX,  see  [7], 


Theorem  1  Let  X  be  an  n  x  m  (n  >  m)  complex  Gaussian 
random  matrix  distributed  asX.  ~  CN(0,  Si  ®  S2),  where 
Si  €  Cnxn  and  S2  £  Cmxm  are  positive  definite  Hermi- 
tian  matrices.  Then  the  density  function  of  S  =  X7/  AX  is 
given  by 


f(S ) 


_ t _ Mpf  c)n-m 

Crro(n)(detSi)™(detS2A)^  ; 

x0F^m\B,-E^S),  (4) 


where  A  £  <C”xn  is  a  positive  definite  Hermitian  matrix 
and  B  =  A-1/2 . 


where  X  £  Cmxm,  Y  £  C"x"  and  n  >  m.  Moreover, 
k  =  (ki, . . .  ,  km)  denotes  a  partition  of  the  integer  k  with 
ki  >  ■  ■  ■  >  km  >  0  and  k  =  ki  +  ■  ■  ■  +  km,  and 
denotes  summation  over  all  partitions  k  of  k.  Tire  complex 
zonal  polynomial,  CK(X),  of  a  complex  matrix  X  defined 
in  [3]  is 

Ck(X)=X[k](1)x[k](X),  (2) 

where  X[k](1)  is  the  dimension  of  the  representation  [k]  of 
the  symmetric  group  given  by 


Tire  distribution  of  the  matrix  S  is  denoted  by 
CQn  m(A,  Si,  S2).  It  should  be  noted  that  the  density  given 
in  Theorem  1  is  different  from  the  one  given  in  [4],  From 
this  generalized  density  we  can  easily  derived  other  well- 
known  densities.  Special  cases  of  density  (4)  are: 

(i)  If  A  =  then  the  density  of  S  =  X7,X  ~ 
i,S2)  is  given  by 

/(S)  =  Crm(„)(detS\)m(detSi)4^>"-" 

(5) 


X[«c](!)  =  k- 


n ?<j(h  -  kj  -  i  +  j) 

HT=1  (h+m-i)\ 


(ii)  If  A  =  In,  Si  =  In  and  S2  =  S,  then  S  =  XHX 
is  said  to  have  a  complex  Wishart  distribution,  denoted  by 


CWm(n,  £),  witli  density 


/(5)  =  Cr„(„xdetE).(det5)"'m^-s~l5>’  (6) 

where  etr  denotes  the  exponential  of  the  trace,  etr(-)  = 

exp(tr(-)). 

(Hi)  If  A  =  In,  £1  =  In  and  £2  =  <r2/m,  then  S  = 
XHX  ~  CWm(n,  cr2Im)  is  a  complex  Wishart  matrix  with 
density 


(2 \-mn  /  i  \ 

m  =  W4-7  (det  5)n_m  etr  5  .  (7) 

cr m(n)  V  &  ) 

The  next  theorem  gives  the  joint  eigenvalue  density  of 
quadratic  forms  on  complex  random  matrices.  This  theorem 
is  one  of  the  key  contributions  of  this  paper.  Moreover,  from 
this  theorem  we  can  easily  derived  other  joint  eigenvalue 
densities  of  complex  random  matrices,  see  [7], 

Theorem  2  Consider  themxm positive  definite  Hermitian 
matrix  S  =  XHAX  ~  CQn>m(A,Y,i,Y,2),  where  A  £ 
Cnxn  is  a  positive  definite  Hermitian  matrix.  Then  the  joint 
density  of  the  eigenvalues,  Ai  >  A2  >  •  •  •  >  Am  >  0.  ofS 
is 


/(A)  = 


nm(m- l)(det£2A)-"  ft  fr,. 

Crm(n)Crm(rn)(det£1)"‘  W  ^  k 


^vCk(5)Cb(-£2-1)Ck(  A) 

k\CK(In)CK(Im) 


At)2 

(8) 


where B  =  A  1/,2£j XA  1/2. 

Special  cases  of  Theorem  2  are: 

(i)  If  A  =  In,  then  the  joint  eigenvalue  density  of  S  = 
XHX  ~  CQnm(In,  Ei,  £2)  is  given  by 


/( A)  = 


7T 


™(m-1)(det  E2)' 


CTm  (n)CTm  (m)  (det  E  i ) 

OO 

*££ 


nAr”n<A*-A')2 


fc=l  k<l 

CK(  E1-1)C'«(-E2-1)C'«(A) 


k=  0  k 


k\CK(In)CK(Im) 


(9) 


(ii)  If  A  =  In,  Ei  =  In  and  E2  =  E,  then  the  joint  eigen¬ 
value  density  of  the  complex  Wishart  matrix  S  =  XHX  ~ 
CWm(n,  E)  is  given  by 


/(A) 


7rm(m-l)(det  £)-n 

crro(n)crro(m) 


k=  1  k<l 


Xo-FJ)m)(— E-1>  A). 


(10) 


(Hi)  If  A  =  In,  Ei  =  In  and  £2  =  a 2Im,  then  the 
joint  eigenvalue  density  of  the  complex  Wishart  matrix  S  = 


XHX  ~  CWm(n,a2Im)  is  given  by 

^mim—  1)  (fr2\  —  nm  771  771 

/(A)  = 


exp(-^EAfc) 


(id 


Next,  we  extend  the  smdy  of  complex  central  Wishart 
distribution  to  the  singular  case,  where  0  <  n  <  m  and 
n,m  £  Z.  Thus  the  rank  of  S  £  Cmxm  is  n.  If  0  <  n  <  m, 
then  the  density  does  not  exist  for  S  ~  CWm(n,  E)  on  the 
space  of  Hermitian  m  x  m  matrices,  because  S  is  singu¬ 
lar  and  of  rank  n  almost  surely.  However,  it  can  be  shown 
that  the  density  does  exist  on  the  (2 mn  —  n2) -dimensional 
manifold,  CSm>n,  of  rank  n  of  positive  semidehnite  m  x  m 
Hermitian  matrices  S  with  n  distinct  positive  eigenvalues. 
Moreover,  the  set  of  all  m  x  n  matrices,  E\,  with  or¬ 
thonormal  columns  is  called  the  Stiefel  manifold,  denoted 
by  CVn>m.  Thus, 

CVn>m  =  {E].(m  x  n );  E?Ei  =  /„}. 


The  application  of  the  next  two  theorems  is  another  key  con¬ 
tribution  of  this  paper.  The  first  of  these  theorems  gives  the 
complex  singular  Wishart  density,  see  [7], 


Theorem  3  The  density  of  S  =  XHX  ~  CWm(n,  E)  on 
the  space  CSm,n  of  rank  n  (0  <  n  <  m)  positive  semidefi- 
nite  m  x  m  Hermitian  matrices  is  given  by 


f(S)  = 


Tn(n—m) 


Cr„(n)(det£)’ 


etr  (-£-1S)  (det  A)r 


(12) 


where  S  =  EiAEf1 ,  E\  £  CVn,m  and  A  = 
diag(Ai , . . .  ,A„). 


The  second  of  these  theorems  gives  the  joint  eigenvalue 
density  of  a  complex  singular  Wishart  matrix,  see  [7], 

Theorem  4  Let  S  =  XHX  ~  CWm(n,  E)  with  0  <  n  < 
m  so  that  S  is  an  m  x  m  positive  semidefinite  Hermitian 
matrix  with  n  positive  eigenvalues.  Then  the  joint  density  of 
the  eigenvalues ,  Ai , . . .  ,  Xn ,  of  S  is 


/(A) 


7r»(n-1)(detS)-” 

Cr„(n)CT„(m) 


IIAr"Il(A*-A;)2 

k= 1  k<l 


XoitfVE-SA), 


(13) 


where  S  =  EiAE^ ,  E\  £  CVn,m,  and  A  = 

diag(Ai , . . .  ,A„). 


In  the  next  section,  we  shall  use  these  densities  to  evalu¬ 
ate  the  two  most  important  information-theoretic  measures, 
namely,  the  ergodic  channel  capacity  and  the  capacity  ver¬ 
sus  outage  of  MIMO  Rayleigh  distributed  channels. 


3.  MULTIPLE-ANTENNA  SYSTEMS 


In  recent  years,  multiple-antenna  techniques  have  become 
a  pervasive  idea  that  promises  extremely  high  spectral  effi¬ 
ciency  for  wireless  communications.  The  basic  information 
theory  result  reported  in  the  pioneering  papers  by  Foschini 
and  Gans  [2]  and  Telatar  [12]  showed  that  enormous  spec¬ 
tral  efficiency  can  be  achieved  through  the  use  of  multiple- 
antenna  systems,  provided  that  the  complex-valued  propa¬ 
gation  coefficients  between  all  pairs  of  transmitter  and  re¬ 
ceiver  antennas  are  statistically  independent  and  known  to 
the  receiver  antenna  array.  However,  in  practical  systems 
the  channel  coefficients  from  two  different  transmitter  an¬ 
tennas  to  a  single  receiver  antenna  can  be  correlated  (at  the 
transmitter  end  with  covariance  matrix  Et)  and/or  from  a 
single  transmitter  antenna  to  two  different  receiver  anten¬ 
nas  can  be  correlated  (at  the  receiver  end  with  covariance 
matrix  £r).  Such  channel  correlation,  which  degrades  ca¬ 
pacity  (see  Chuah  et  al.  [1]),  depends  on  the  physical  pa¬ 
rameters  of  the  MIMO  system  and  the  scatterer  charac¬ 
teristics.  The  physical  parameters  include  the  antenna  ar¬ 
rangement  and  spacing,  the  angle  spread,  the  angle  of  ar¬ 
rival,  etc.  One  of  the  objectives  of  this  paper  is  to  evaluate 
this  capacity  degradation  for  the  correlated  channel  matrix 
H  CN (0,  Er  <g>  £*)  with  nr  >  nt,nr  =  nt  and  nt  >  nr. 

The  complex  signal  received  at  the  jth  output  can  be 
written  as 

nt 

Vi  =  ^2hUxi  +vi’  04) 

i=  1 

where  h,j  is  the  complex  channel  coefficient  between  input 
i  and  output  j,  Xt  is  the  complex  signal  at  the  zth  input  and 
Vj  is  complex  Gaussian  noise  with  unit  variance,  as  shown 
in  Figure  1.  The  signal  vector  received  at  the  output  can  be 
written  as 


i.e.,  in  vector  notation, 

y  =  Hx  +  v,  (15) 

where  y,v  E  Crtr ,  H  £  <C”rXn*,  x  £  C”‘  and  v  ~ 
CN( 0,  In,.)-  The  total  power  of  the  input  is  constrained  to 
P , 

£{xHx}  <  p  or  tr  £{xxH}  <  p. 

We  shall  deal  exclusively  with  the  linear  model  (15)  and 
derive  the  capacity  of  MIMO  channel  models  in  this  section. 


Fig.  1.  A  MIMO  communication  system. 


3.1.  Ergodic  channel  capacity 

We  assume  that  H  is  a  complex  Gaussian  random  matrix 
whose  realization  is  known  to  the  receiver,  or  equivalently, 
the  channel  output  consists  of  the  pair  (y,  H).  Note  that  the 
transmitter  does  not  know  the  channel  and  the  input  power 
is  distributed  equally  over  all  transmitting  antennas.  More¬ 
over,  if  we  assume  a  block-fading  model  and  coding  over 
many  independent  fading  intervals,  then  the  Shannon  or  er¬ 
godic  capacity  of  the  random  MIMO  channel  [12]  is  given 
by 

C  =  £h  {logdet  [Int  +  (p/nt)HHH)} 

=  £h  {logdet  (J„r  +  (p/n*)HH")},  (16) 

where  the  expectation  is  evaluated  using  a  complex  Gaus¬ 
sian  density,  i.e.,  H  ~  CiV(0,£r  ®  £*).  Let  S  =  HhH. 
Then  the  channel  capacity  can  be  written  as 

C  =  £s  {logdet  (4,  +  (p/nt) S)}  ,  (17) 

where  the  expectation  is  evaluated  using  the  density  given  in 
(5).  Note  that  if  the  channel  is  correlated  only  at  the  trans¬ 
mitter  end  and  nt  >  nr  then  S  =  HhH  is  a  complex  sin¬ 
gular  Wishart  matrix  and  its  density  is  given  in  Theorem  3. 
Similarly,  if  the  channel  is  correlated  only  at  the  receiver 
end  and  nr  >  nt  then  S  =  HH/f  is  a  complex  singular 
Wishart  matrix.  For  the  other  situations,  S  is  full  rank  and 
the  distributions  are  chosen  according  to  the  spatial  correla¬ 
tion  types,  e.g.,  (5)—  (7)  (see  Table  1  in  [7]). 

Let  Ai  >  •  •  •  >  \nt  be  the  eigenvalues  of  S  and  A  = 
diag(  Ai , . . .  ,  \nt )  •  Then  the  capacity  can  also  be  computed 
using  the  joint  eigenvalue  density,  /(A),  i.e., 


Again,  the  expectation  in  ( 1 8)  is  evaluated  using  the  density 
according  to  the  spatial  correlation  types,  e.g.,  (9)  —  (1 1)  and 
(13). 


3.2.  Capacity  versus  outage 

Another  useful  information- theoretic  measure  is  the  ‘x  per¬ 
cent  outage’  which  is  defined  to  be  the  minimum  mutual 
information  that  occurs  in  all  but  x  percent  of  the  channel 
instantiations.  In  other  words,  if  we  measure  the  mutual 
information  of  the  channel  many  times  —  in  many  instanti¬ 
ations  of  the  random  channel  —  we  would  find  that  a  mutual 
information  greater  than  5%  outage  would  occur  95%  of  the 
time  (x  =  5). 

Let  us  formulate  the  outage  capacity.  For  a  given  instan¬ 
tiation  of  H,  if  the  receiver  knows  the  channel,  the  mutual 
information  X(x;  (y,  H))  is  given  by 


Here  Xout  =  out(P0ut)  means  that,  in  any  instantiation  of 
H  from  the  ensemble,  we  shall  obtain  a  mutual  information 
X  less  than  Xout  with  probability  Pout-  Using  the  density  of 
the  mutual  information  we  can  give  another  expression  for 
the  ergodic  capacity  as  follows: 

POO 

C=  Xq(X)dX.  (23) 

Jo 

In  other  words,  the  mean  of  the  mutual  information  gives 
the  ergodic  capacity. 

4.  COMPUTATION  OF  THE  CAPACITIES 


Z(x;(y,H)) 


log  det 
nt  / 

k= 1  ' 


Int  +  ~  KH  H 


rh 


1  +  -A* 
nt 


(19) 


where  Ai, . . .  ,  \Ut  are  the  eigenvalues  of  S  =  HffH.  The 
density  q(X)  of  the  mutual  information  X  in  equation  (T9) 
over  the  ensemble  of  instantiations  of  H  can  then  be  written 
as 


where  6  is  the  Dirac  delta  function,  see  Simon  and  Mous- 
takas  [11].  The  following  theorem  gives  this  density  in 
terms  of  the  density  of  a  single  unordered  eigenvalue  of 

S  =  HhH. 


Theorem  5  Let  H  be  an  nr  x  random  matrix  channel 
and  g{  A)  be  the  density  of  a  single  unordered  eigenvalue  of 
S  =  H^H  Then  the  density  of  the  mutual  information, 
q(X),  is  given  by 


where  p  is  the  total  transmitter  power. 


We  can  also  express  the  distribution  Q(X)  of  the  mutual 
information  X  as  follows: 


q(X)  dX. 


(22) 


The  unique  inverse  function  of  Q(X)  is  called  the  out¬ 
age  mutual  information  and  is  denoted  by  out(P0ut)  where 
Pout  =  <9 (lout).  In  other  words,  we  define  out(_Pout)  such 
that 

Pout  =  Out  (Pout) 


hi  this  section,  we  evaluate  the  capacities  for  both  correlated 
and  uncorrelated  Rayleigh  2x2  fading  channels  ( nr  = 
nt  =  2).  First,  we  consider  a  channel  with  correlation  at 
both  transmitter  and  receiver  ends.  A  typical  example  of 
this  situation  is  a  communication  between  two  mobile  lap¬ 
tops  in  poor  scattering  conditions.  In  addition,  due  to  phys¬ 
ical  size  constraints  it  is  more  difficult  to  space  the  antennas 
far  apart  at  the  laptops.  These  constraints  lead  to  channel 
correlation  at  both  transmitter  and  receiver  ends,  i.e.,  H  ~ 
CN( 0,£r  <g>  St)  and  S  =  ET^H  ~  CQnr.,nt  (Inr,  st)- 
The  ergodic  capacity  is  given  by 


(det  St) ' 


(detSr)2 

[l  +  2^2])  ^ 


J™  J ^  dX2  d\i  (log  [l  +  |Ai] 


+ 


A  2f 


y-  y^  X[k](^v  *)  X[«](-^t  1)  X[k](A) 

^4'(fc1-fc2  +  i)r(fc1  +  2)r(fe  +  i)‘ 


Second,  we  consider  a  channel  with  correlation  only  at  the 
transmitter  end.  A  typical  example  of  this  situation  is  an 
uplink  communication  from  a  mobile  unit  to  a  base  station. 
Here,  the  antennas  at  the  base  station  can  be  spaced  suffi¬ 
ciently  far  apart  to  achieve  uncorrelation  at  the  receiver  end 
but,  due  to  physical  size  constraints,  it  is  more  difficult  to 
space  file  antennas  far  apart  at  the  mobile  unit  (transmitter 
end),  which  leads  to  correlation  at  the  transmitter  end,  i.e., 
H  ~  CN( 0,  In„  ®  St)  and  S  =  HhH  ~  CWnt  ( nr ,  St)  is 
a  complex  Wishart  matrix.  It  should  be  noted  that  the  joint 
eigenvalue  density  of  a  Wishart  matrix  depends  on  the  pop¬ 
ulation  covariance  matrix  £*  only  through  its  eigenvalues, 
Vi, . . .  ,  vnt ,  i.e., 

0P0(nt)  (-£t-\  A)  =  oPo(nt)  (-T-1,  A) , 

where  T  =  diag(rti,...  ,vnt).  Let  nt  =  2  and  T-1  = 
diag(oi,  02 ).  Then  from  [6]  we  have 


oP0(2)(-T-1,A)  = 


(02  —  ai)(Ai  —  A2) 


X  |"g— 01A1— 0,2X2  g— (X1A2— 0,2X1' 


(24) 


where 


Pout  —  Q(2out). 


The  ergodic  capacity  is  given  by 
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transmitter  and  receiver  ends,  respectively. 
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In  this  paper,  we 

•  study  the  densities  of  quadratic  forms  on  complex  random  matrices 

•  study  the  eigenvalue  densities  of  quadratic  forms  on  complex  random  matrices 

•  show  that  these  densities  are  represented  by  complex  hypergeometric  functions 
of  matrix  arguments 

•  show  that  the  above  can  be  expressed  in  terms  of  complex  zonal  polynomials 

•  derive  a  close  form  ergo  die  capacity  and  capacity  versus  outage  formula  for 
spatially  correlated  MIMO  channels 

•  show  how  the  channel  correlation  degrades  the  capacity  of  the  communication 
system 
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Figure  1:  A  MIMO  communication  system. 

The  signal  vector  received  at  the  output  can  be  written  as 

y  =  Hx  +  v ,  where  y,v  G  Cn\  x  e  Cnt ,  H  e  CUrXnt.  (1) 

The  total  power  of  the  input  is  constrained  to  p, 

£{xhx}  <  p  or  ti£{xxH}  <  p. 
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The  ergodic  capacity  of  MIMO  channel  is  given  by 

C  =  max£H{^(x;y)} 

/(*) 

=  £h  {log det  (Int  +  (p/m) HWH) } 

=  £s  {log det  (Int  +  (p/nt) S)} 

=  £a  /  log  [1  +  (p/nt)\k\ 

JH 

=  8 **  {*°S(1  +  (P/Ut)^k)}  ,  (2) 

k= 1 

where  the  expectation  is  evaluated  using  the  density  according  to  the  spatial 
correlation  types,  see  the  table  in  the  next  slide. 

If  S  =  HhH  ~  CQnrnt(Inr,  Er,  Et)  then  the  channel  is  spatially  correlated 
Rayleigh  distributed,  where  Er  and  Et  represent  the  channel  correlation  at  the 
receiver  and  transmitter  ends,  respectively. 
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Spatial  correlation 

Channel  distributions 

types 

nr  >  nt 

s 

*3 

II 

'S 

C+- 

nt  >  nr 

Uncorrelated 

(i.i.d) 

H  ~  CN(0,  Inr  ®  <t24,)= 
S  =  H"H 

S  ~  CW„,(nnA,) 

<= 

H  ~  CN( 0,  cr2/„r  ®  Int) 

S  =  HH" 

=  S  ~CWnr(nt,a2Inr) 

Correlated  at 

the  transmitter 

H  ~  CN(0,  IHr  ®  Et) 

g  _ 

S  ~  CWnt(nr,S() 

H  ~  CN(0,  Inr  ®  Et) 

g _ 

S~CWnt(nnEt) 

Correlated  at 

the  receiver 

H  ~  CN(0,  Er  ®  /„,) 

g  _ 

S  ~  CWnr(nt,  Er) 

H  ~  CN( 0,  Er  ®  /„,) 

g  _ 

S  ~  CWnr(nt,  Er) 

Correlated  at  both 

transmitter  and  receiver 

H~CiV(0,Er®Et)  = 

g  _  jjff 

S  ~  CQnr  nt(I„r,  Er,  Et) 

<= 

H  ~  CN(0,  Er  ®  Et) 

g  _ 

CJ  /T  /”»  /  T  V'"l  Y"1  \ 

~  CC^nrjnt(int,  2^r,  2^J 
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Signal-to-Noise  Ratio  (in  dB) 


Figure  2:  Ergodic  capacity  vs  SNR  for  nt  =  2  and  nr  =  2,  i.e.,  Rayleigh  2x2  channel 
matrices.  Note  that  Cu,  Cc  and  Ccc  denote  the  capacity  of  the  three  environments: 
uncorrelated,  correlated  at  the  transmitter  end  and  correlated  at  both  transmitter 
and  receiver  ends,  respectively. 
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